Triplet superconductivity in a ferromagnetic vortex 
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We argue that odd-frequency triplet superconductivity can be conveniently realized in hybrid 
superconductor-ferromagnet (SF) structures with a ferromagnetic vortex. We demonstrate that 
due to proximity-induced long-range triplet pairing such SFS junctions can sustain appreciable 
supercurrent which can be directly measured in experiments. 



A normal metal (N) sandwiched between two super- 
conductors (S) can become superconducting as a result 
of penetration of Cooper pairs from the superconduct- 
ing electrodes. The range of penetration is set by the 
so-called thermal length which can easily reach sev- 
eral micrometers at sufhciently low temperatures [iHll- 
The situation changes drastically if the normal metal is 
replaced by a ferromagnet (F) . The quantum mechanical 
exchange interaction on the F-side then destroys conven- 
tional spin-singlet Cooper pairs within a few nanometers 
(the so-called paramagnetic effect) Experiments to 
determine actual superconducting penetration depths in 
ferromagnets intensified more than a decade ago, when 
techniques were developed to fabricate hybrid nanoscale 
SF structures with well controlled geometries. Several 
groups 0-0] reported an unexpectedly strong influence 
of superconductors that stimulated new theoretical ef- 
forts in a search for a sustainable superconductivity that 
is compatible with the exchange interaction. During 
the last decade several theoretical mechanisms were sug- 
gested [iMl, some of which were successfully realized 




FIG. 1: (Color online) SFS junction formed by two supercon- 
ducting electrodes connected via ferromagnetic vortex. 



structures have attracted a lot of attention due to their 



experimentally 17h22||. A recent comprehensive review 
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of the status of the field was given in Ref. 

Common to all mechanisms of long range proximity 
effect in ferromagnets is the generation of triplet super- 
conductivity within highly inhomogeneous ferromagnetic 
regions adjacent to superconductors. The systems stud- 
ied up to date include intrinsically inhomogeneous ferro- 
magnets , half- metallic ferromagnets with spin-active 
FS interfaces [Ulii], and engineered multilayers consist- 



ing of magnetic and non-magnetic materials 

In this letter we address a different situation of 
proximity-induced long range triplet pairing in ferromag- 
nets with magnetic vortex structure. Magnetic vortices 
are stable in systems intermediate between very small, 
10 nm scale magnets, which behave as single giant spins, 
and macroscopic magnets with dimensions exceeding ~ 1 
yum. The magnetic structure in such mesoscopic mag- 
nets is the result of a competition between exchange, 
anisotropy, and dipolar energies and depends strongly on 
their shape. The latter property allows magnetic nano- 
engineering using modern nanolithography [24J opening 
possibility of investigating SFS structures with differ- 
ing magnetic structures. Recently mesoscopic magnetic 



remarkable transport properties |25l - l27l |. as well as the 
prospect of technological applications for magnetic stor- 
age of information of unprecedented density [28] . Below 
we will demonstrate that mesoscopic ferromagnetic struc- 
tures can turn superconducting if attached to supercon- 
ducting electrodes. 

The model and quasiclassical formalism. We will con- 
sider a ferromagnetic film of thickness d located in the 
xy plane with magnetization forming a vortex. This film 
is partially covered by two superconducting electrodes 
thus forming an SFS contact as it is shown in Fig. [T] 
Our main goal is to analyze superconducting correlations 
that penetrate into a ferromagnetic vortex from the elec- 
trodes. In order to accomplish this goal we will employ 
the quasiclassical Usadel equations 0, for energy- 
integrated matrix Matsubara- Green functions G. E.g. 
in the ferromagnet with diffusion constant D these equa- 
tions read 
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are 4x4 matrices in Nambu and spin spaces. Their 
commutator in Eq. ^ and below is denoted by square 
brackets. Accordingly, G, F, F+ and (5+ are 2x2 matri- 
ces in the spin space, a;„ = TrT{2n + 1) is the Matsubara 
frequency, h is the exchange field in the ferromagnet and 
(T — ((7i, i5'2, (73) represents the Pauli matrices in the spin 
space. The same equations ([1]), ([2]) hold also for super- 
conducting electrodes, one should only replace D by the 
diffusion constant in the corresponding electrode. The 
superconducting order parameter A equals to zero in the 
ferromagnet, while in two superconducting terminals it is 
respectively A = Aiexp(zx/2) and A = A2exp(— zx/2) 
with real Ai_2 and x being the superconducting phase 
difference across our SFS junction. 

Equations ([T]) should be supplemented by appropri- 
ate boundary conditions at each of the two SF-interfaces 
which account for electron transfer across these inter- 
faces. In what follows we will assume that there exist 
tunnel barriers at both SF interfaces with the correspond- 
ing tunneling resistances ri^2- In the tunneling limit it 
suffices to employ Kuprianov-Lukichev boundary condi- 
tions fs^l at each SF- interface. E.g., at the interface be- 
tween the first superconducting electrode {z > 0) and the 
ferromagnet {z < 0) these boundary conditions read 

2naGFd,GF^[GF,GsA, (3) 

where Gp and Gs-^ are respectively the Green functions 
at the F- and S-sides of the first interface and a is the 
Drude conductivity of a ferromagnet. Analogous bound- 
ary conditions hold for the second SF-interface. 

Long-range triplet pairing in a ferromagnetic vortex. 
The presence of tunnel barriers at both SF-interfaces ef- 
fectively implies weak electron tunneling regime in which 
case the proximity effect remains small and it suffices to 
linearize Usadel equations in the ferromagnet as 

DV'^F-2ujnF~i{F,h{r)&} = 0. (4) 

In Eq. (HI) we restrict Matsubara frequencies to be pos- 
itive a;„ > and denoted the anticommutator by curly 
brackets. A similar equation holds for the function F'^. 

In general magnetization patterns in thin ferromag- 
netic films depend on the film geometry and are influ- 
enced by the following trade-off. On one hand, magneto- 
static energy minimum is reached provided the film mag- 
netization remains in-plane. On the other hand, in some 
regions, such as, e.g., vortex cores, local magnetization 
can go out-of-plane in order to minimize the exchange en- 
ergy. As the magnetic core radius typically remains small 
as compared to the superconducting coherence length, 
in the following we will assume that magnetization lies 
in-plane everywhere in the ferromagnet, see Fig. [TJ In 
sufficiently thin films the exchange field h depends only 
on in-plane coordinates {x,y) and can be represented as 
h = (/icos6', /isin6') where 9 = 6{x,y). In this case the 
spin structure of the anomalous Green function F inside 



the ferromagnet can be chosen in the following form 

F^Fo + &mFh + &[e„m]Ft, (5) 

where Fq describes the singlet pairing component, while 
Fh and Ft correspond to two different triplet components. 
In Eq. ([5]) we also introduced in-plane and normal to the 
plane unity vectors m = h/h and e^. Combining Eqs. 
^ and (jll) we arrive at the following equations for the 
above components: 

DV^Fo - 2ajnFo = 2ihFh, (6) 
VF,, = DFtVle + 2D{VpFt, V pff) + 2ihF^, (7) 
VFt = ^DFhVle - 2D{VpFh, VpO), (8) 

where we defined the differential operator 

V = DV^ - D{Vpef - 2ujn (9) 

and distinguished V and Vp as respectively 3d and 2d 
(in-plane) gradient operators. 

Note that Eqs. ([6]) and (O contain the exchange field 
h thus providing the characteristic length scale both for 
Fq and Fh of order ^ \/D/h. At the same time, 
Eq. ^ does not contain the h-teim and, hence, typi- 
cal variations of Ft occur on a much longer length scale 
~ \/D/T :s> ^h- This observation illustrates the dif- 
ference between the two triplet components Fh and Ft 
and constitutes the essence of the long range proximity 
effect in SFS structures: while the components Fq and 
Fh decay already in the vicinity of an SF-interface, the 
triplet component Ft survives deep inside the ferromag- 
net provided the temperature remains sufficiently low. 

Before turning to the solution of Eqs. ([B])-® let us 
perform some further simplifications. Firstly, we will ne- 
glect both magnetic anisotropy and stray field effects. In 
this case outside the magnetic vortex core the function 9 
obeys the equation 

= 0, (10) 

which allows to drop the first terms in the right-hand 
side of Eqs. ([7]) and ([8|). Secondly, we will assume the 
ferromagnetic film to be sufficiently thin d , in which 
case the dependence of the long-range triplet component 
Ft on the coordinate z can be neglected. Then, integrat- 
ing Eq. ([8|) over z we obtain 

VpFt^-2D{VpFh,Vp9), Fh = -j ^Fhdz, (11) 

where Fh is the average value of Fh component over the 
ferromagnetic film thickness and Vp is defined by Eq. (|9]) 
with ^ V^. 

Eq. (fTT|) accounts for diffusion of the long-range triplet 
component Ft across the ferromagnet with nonuniform 
in-plane magnetization. It demonstrates that non-zero 
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Ft is generated in the parts of the ferromagnet where 
both Vpd and V pFh differ from zero. The condition 
Vp^ 7^ obviously holds everywhere in the ferromag- 
netic plane since the magnetization remains non-uniform 
there. As for the averaged component Fh: it vanishes 
together with its gradient at distances exceeding S^h 
from SF- interfaces. In the immediate vicinity of such 
interfaces Fh is non-zero, but its in-plane gradient re- 
mains small because in the main approximation it only 
depends on the absolute value of the exchange field h, 
cf. Eq. d?]). The gradient V pFh becomes appreciable 
only in the region of the ferromagnet just below the edge 
of the superconducting film where Fh changes abruptly. 
With this in mind we arrive at the following result for 
the long-range triplet component 



had Tfe 

fc=l,2 



(12) 

which holds inside the ferromagnetic film. Here Fs^ is 
anomalous Green function in the bulk of the fc-th super- 
conductor and riij, is the outer unity vector normal to 
the superconducting plane Sk (see Fig. [ij. Integration 
contours Ik in Eq. (|12|) are lines in the xy plane corre- 
sponding to the edge of the superconductor Sk and in the 
ferromagnet kernel P^'^ obeys the equation 



(13) 



with boundary conditions dP^'^ jdn — 0. We also note 
that Eq. ([T2|) can easily be generalized to the case of 
arbitrary Vp6' not obeying Eq. pO)) . 

Triplet pairing and Josephson effect. As triplet pairing 
amplitude can survive deep in the ferromagnet, at suffi- 
ciently low temperatures our SFS junction can sustain 
appreciable supercurrent which is converted from singlet 
to triplet and back in the vicinity of SF-interfaces. In 
order to evaluate this supercurrent we will employ the 
standard expression for the current density 

j = Im Sp[FVF+ ~ F+\7F], (14) 

w„>0 



2e 



where the trace is taken over the spin degree of freedom. 
Combining Eqs. ([S|), (IT^ with (|14l) we recover the sinu- 
soidal current-phase relation I{x) = Ic sinx with 

^ (15) 
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Note that in the course of our derivation we always as- 
sumed the proximity effect to be sufficiently weak. This 
assumption is satisfied under the condition 
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(16) 
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FIG. 2: (Color online) Typical spatial distribution of the long- 
range superconducting triplet component Ft induced in the 
ferromagnetic disk with vortex-like magnetization by one su- 
perconducting electrode {x > R/2, z > Q) with real A. 



Eq. ([15]) together with its validity condition repre- 
sents the central result of our analysis which fully deter- 
mines the Josephson critical current of an SFS junction 
with a ferromagnetic vortex. Actually this result applies 
not only to vortex configurations but also to a broader 
class of non-uniform magnetization patterns. 

Let us now assume that our ferromagnetic film has the 
form of a disk with radius R and vortex-like magnetiza- 
tion pattern with vortex core located in the disc center. 
Then the function equals to tp + 11/2 for clockwise or 
</j — 7r/2 for counterclockwise magnetization, where (p is 
the azimuthal angle (see Fig. [T]). Eq. ([TU|) is fulfilled 
in this case. Remarkably, Eq. P3|) yields exactly the 
same result for quite different magnetization patterns: 
vortex-like [0 = ip ± 7r/2), antivortex-like [9 ~ —p or 
Q — —if -I- tt) and hedgehog-like {9 — p or 9 = p + tt) 
states. This property holds since the function Vp6' re- 
mains the same (up to a sign) for all these magnetization 
patterns. Note, however, that for the last two patterns 
stray magnetic field is not confined to the disc center and 
may influence superconductivity in the electrodes. 

For illustration, typical spatial profile of the long-range 
superconducting triplet component Ft induced by one su- 
perconducting electrode in the ferromagnet with a vortex 
is schematically depicted in Fig. [2] As it was expected. 
Ft is most efficiently generated close to the edge of a su- 
perconductor where the scalar product |(Vp6'(p), n/(p))| 
reaches its maximum values. Provided the proximity ef- 
fect remains weak, the total value of Ft is given by a 
superposition of independent contributions from two su- 
perconducting electrodes, cf. also Eq. (IT2t . 

As one can observe in Fig. [21 the long-range triplet 
component Ft penetrating into the ferromagnet can take 
both positive and negative values. Thus, depending on 
the magnetization pattern inside the ferromagnetic film 
it is possible to realize both zero- and 7r-junction states 
in our structure. The latter regime can be reached, e.g., 
by implementing certain asymmetry in SF contacts. 

We further consider a symmetric situation, set Ai_2 = 
|A| and assume that the relevant Thouless energy exh ~ 
D/{2R)'^ remains smaller than the superconducting gap 
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In summary, we demonstrated that long-range triplet 
superconductivity can coexist with a ferromagnetic vor- 
tex and evaluated the supercurrent across SFS junctions 
containing such vortex. For properly chosen system pa- 
rameters the effect is well in the measurable range and 
can be directly tested in future experiments. This work 
was supported in part by DFG, by RFBR under grant 
09-02-00886 and by British EPSRC grant EP/F01689/1. 



FIG. 3: (Color online) /c versus temperature (normalized by 
the critical temperature Tc) in SFS junctions containing a 
ferromagnetic vortex at different values of R. The edges of 
superconducting electrodes (contours ^1,2) are chosen to coin- 
cide with straight lines y = ±i?/2. 



|A|. Then in the limit T < exh from Eq. ((T5l) we find 

Ic - D^ETy^/iedh^^r^a), (17) 

while at intermediate temperatures exh ^ T ^ |A| the 
Josephson current follows the standard exponential de- 
pendence on temperature 
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-exp (^~L^/2^TT/DJ 



(18) 



where L is an effective distance between the two SF 
contacts which depends on geometry details (obviously 
L = 2R for small area contacts). For illustration the 
Josephson critical current Ic is also plotted in Fig. |3] as 
a function of temperature for different values of R. 

Our result for Ic in SFS systems turns out to be by 
the factor ~ e'^^^/h'^ smaller than that for conventional 
diffusive SNS junctions with identical geometry, cf., e.g., 
d, [s^. The critical current of our SFS structure can 
further be increased by a proper choice of the system 
parameters. For a simple estimate of possible maximum 
values of Ic let us employ Eq. (IT7|) at the border of its 
applicability range (ITS]) . Then for T ^ exh and d > 
we obtain 



D^a/{eR^dh) ^ (a/d)2exh/(ei?Ar), 



(19) 



where Rn is the normal state resistance of the ferromag- 
netic film between two superconducting electrodes. This 
estimate is also supported by our independent calcula- 
tion (not presented here) which yields contributions to 
Ic on 1/h in higher orders in barrier transmissions. Eq. 
(fT9| demonstrates that for d^ one can expect to reach 
values of Ic only by the factor ~ smaller that the 

absolute maximum Ic £xh/ei?Ar achieved for SNS junc- 
tions [sij. Actually, the latter maximum value can also 
be reached, but only for extremely thin films d ^ (cf. 
Eqs. dnl), dH])) with large values oi Rn- 
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